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(57) Abstract: A method for estimating a receiver's location (X) in a wireless communication environment (RN) having several 
channels. Each channel has at least one signal parameter (V) that varies with location (X) differently from the other channels. A 
set of calibration data (CD) is determined for each calibration point, each set comprising the location (X) and at least one measured 
signal parameter (V) for each of several channels. The calibration data (CD) serve as a basis for a statistical model (SM) of the signal 
parameters (V) versus a receiver's location. A set of observed signal parameters (CO) is determined, the set comprising at least one 
signal parameter (V) for each of several channels at the receiver's location (X). A location estimate (LE) approximating the location 
(X) of the receiver (R) is determined on the basis of the statistical model (SM) and the set of observed signal parameters (CO). 
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Location estimation in wireless telecommunication networicis 

Background of the invention 

The invention relates to metiiods and equipment for estimating a 
receiver's location in a wireless telecommunication environment, ie one or 

5 more networks which may be radio, microwave or optical networks. The one or 
more networks communicate at a plurality of channels simultaneously. Such a 
location estimation can be used to provide a wide variety of location- 
dependent services. 

US patent 6 1 12095 to Mati Wax et al. discloses a method for pro- 

10 viding a set of likely locations of a transmitter in a cellular network, such as 
AMPS or CDMA. A problem with the technique disclosed in the Wax patent is 
that it requires additional hardware at the network side, such as an antenna ar- 
ray which is equipped to measure an angular direction relative to a base sta- 
tion. In other words, to detennine a mobile station's location, infonmation on 

15 the network infrastructure must be available and the mobile station must 
transmit something to have its location estimated. 

Disclosure of the invention 

An object of the Invention is to solve the above problems. In other 
words, the mechanism according to the invention should be able to estimate a 
20 receiver's location in a wireless telecommunication network even without prior 
knowledge of the network infrastructure (such as the locations of the base sta- 
tions). 

This object is achieved with a method and equipment which are 
characterized by what is disclosed in the attached independent claims. Pre- 
25 fen-ed embodiments of the invention are disclosed in the attached dependent 
claims. 

The invention is based on the surprising idea that it is possible to 
estimate a receiver's location with reasonable confidence without knowledge 
of the infrastnjcture of the receiver's wireless environment, ie the network(s) 

30 received by the receiver. For example, the technique disclosed in the above- 
referenced Wax patent relies on the cellular network's base station configura- 
tion, including the location of the base stations, it is Indeed surprising that the 
technique according to the invention is feasible. The fact that it is surprising is 
apparent as soon as one walks around with a mobile phone having a field 

35 strength indicator. In some places, a shift of 20 to 30 cm changes the field 
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strength dramatically. Evidently, there must be a vast number of locations with 
near-identical field strength. One would expect that calibrating a location esti- 
mation system requires field strength (or other signal parameter) measure- 
ments at locations very dose to each other, and that huge databases would be 

5 required to store these measurements. Atmospheric conditions, cityscapes 
and network configurations change continuously. At first sight, it would seem 
that the databases will deteriorate rapidly, unless constantly updated. How- 
ever, computer simulations show that a technique based on measurements at 
several channels (frequencies) is surprisingly robust Also, calibration data can 

10 be collected automatically at various conditions. 

One aspect of the invention is a method for estimating a location of 
a receiver In a wireless telecommunication environment, the telecommunica- 
tion ehvlronment comprising a plurality of channels for sirtiultaneous commu- 
nication, each channel having at least one signal parameter that varies with lo- 

15 cation differently from the other channels. The method can be Implemented by 
the following steps: 

1) for each of a plurality of calibration points In the wireless tele- 
communication environment, detemiining a set of calibration data, each set of 
calibration data comprising the location of the respective calibration point and 

20 at least one measured signal parameter for each of several channels at that 
calibration point; 

2) on the basis of the sets of calibration data, maintaining a statisti- 
cal model of the signal parameters of the several channels versus a receiver's 
location in the wireless telecommunication network; 

25 3) measuring at least one signal parameter for each of several 

channels at the receiver; and 

4) estimating the location of the receiver on the basis of the statisti- 
cal model and the measured signal parameters of the several channels at the 
receiver. 

30 Another aspect of the invention is an anrangement for carrying out 

the above method. The arrangement can be embodied as a receiver compris- 
ing means for detenmining sets of observed signal parameters, each set com- 
prising at least one observed signal parameter for each of several channels at 
the location of the receiver. The receiver may itself comprise a location calcu- 

35 iation module for detemiining a location estimate approximating ttie location of 
tiie receiver on the basis of said sets and a statistical model of tiie signal pa- 
rameters of the several channels versus a receiver's location in a wireless 
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telecommunicait|on environment Alternatively, the receiver may convey the 
sets to an external location calculation module. 

The temn 'receiver* means that the device whose location Is being 
estimated does not have to transmit when its location Is being estimated. In 

5 other words, it suffices that the device is making observations of its wireless 
environment. For example, a GSM phone does not have to receive a trafRc 
channel. Rather it maizes observations at ail available frequencies. The device 
may also have, and typically has, transmitting ability, but it is not necessary for 
all embodiments of the invention, and the invention can be used to estimate 

10 the location of a pager or a broadcast receiver. Because transmission capabil- 
ity is not essential to location estimation according to the invention, tiie re- 
ceiver may exploit signal parameters of networks it is not attached to. For ex- 
ampie, a GSM phone attached to one GSM network may exploit signal 
strengtii values of other GSM networks. 

15 The tenn 'environment* means that ttie receiver can receive (make 

observations of) at least one network, but it can receive more tiian one. For 
example, a GSM phone may observe several operators' GSM networks. A 
more advanced receiver may observe many types of networks, such as cellu- 
lar networks and broadcast networks. 

20 A 'wireless' environment means that the one or more networks may 

be radio, microwave or optical networks. Also, the set of networks received by 
the receiver must communicate at a plurality of channels simultaneously, and 
the plurality of ctianneis must comprise a subset of channels such that each 
channel In the subset has at least one signal parameter that varies witti loca- 

25 tion differently from the otiier channels in the subset. This means tiiat several 
channels having signal parameters with near-identical dependence from loca- 
tion, such as channels from a common transmitting antenna, do not nonnally 
give sufficient infonnation for reliable location estimation. Nonnally, signals 
from at least ttiree transmitting stations are required. Examples of suitable 

30 networks are cellular networks (such as GSM, GPRS, UMTS, etc.), broadcast 
networks (analogue audio, DAB or DVB), wireless local-area networks (WLAN) 
or short-range microwave networks, such as Bluetooth. 

A 'location' may have one to three dimensions. A one-dimensional 
presentation of location may be sufficient in trains and the like. Two- or three- 

35 dimensional presentations of location are much more useful, however. In a 
two-dimensional presentation, ttie receiver is assumed to be substantially at 
earth level. Actually tiie height does not matter as long as the calibration data 
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is measured at the same height (siich ias ground level, 13th floor, etc.) as the 
actual observations. Additionally, flie calibration data may comprise a presen- 
tation of time. This means that the wireless environment, le its signal parame- 
ters, vary with time. In other words, the calibration data comprises, in addition 

5 to the signal parameters, one to three location coordinates and, optionally, a 
presentation of time. 

The tenrt 'calibration data', as used herein, comprises calibration 
measurements (le, measured signal values) and the location (and, optionally, 
time) at which the measurements were made. 

10 The temi 'statistical model' means that tiie individual sets of calibra- 

tion data are not needed to calculate ah individual receiver's location. The dif- 
ference between a statistical model and the sets of calibration data can be il- 
lustrated by the following example. Assume tiiat we have a number of {x, y} 
pairs such that there is some dependence between x and y. The y value at a 

15 location X can be calculated on the basis of all the {x, y} pairs. A much faster 
way to predict the value of y given a value of x is to calculate a mathematical 
function y = f(x). In this example, tiie function f Is the statistical model. In otiier 
words, the value of y given a value of x is calculated without reference to tiie 
individual {x, y} pairs. Location estimation on the basis of the statistical model 

20 is faster and requires less storage space tiian location estimation on the basis 
of the individual sets of calibration data. 

The statistical model can have a large variety of different implemen- 
tations, such as probabilistic models, neural networks,, fuzzy-logic systems, 
kernel estimators, support vector machines, decision trees, regression trees, 

25 Kalman filters and other stati'stical filtering methods, wavelets, splines, Induc- 
tive logic programming methods, finite mixture models, hidden Markov models, 
etc. As used in this context, ttie terni 'statistical model' may also refer to a mix- 
ture of several statistical (sub)models. 

The term 'channel' should have a wide interpretation, meaning 

30 more or less the same as a frequency or frequency band. The receiver does 
not have to communicate on the channel, as long as the receiver (or an at- 
tached measuring apparatus) can measure at least one signal parameter of 
that channel. In TDMA systems, each frequency has several timeslots, each of 
which carries one channel. As far as the invention is concerned, all timeslots 

35 having the same frequency give identical information, and any one of them 
can be used as a 'channel'. If the measured signal parameter is signal 
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strength, the receiver does not even have to be able to interpret the contents 
ofthechannei. 

An illustrative but non-exhaustive list of the signal parameters vary- 
ing with location comprises signal strength, timing advance and enior ratio. The 

5 list may ajso comprise the availability of certain channels, but this can be seen 
as a special case in which the signal strength and/or error ratio is quantified to 
a yes/no question. If directional antennas are used, the direction of the radio 
beam(s) can be used as well, thus the measured signal parameters do not 
have to conrespond to a certain channel, but they can be derived values. For 

10 example, a measured parameter set may be or comprise a vector V = |yi , V2, 
V3...] In which V1, V2 etc. are the indices of the best, second best, etc. avail- 
able channel. For the purposes of clarity, however, we will use examples in 
which the signal parameters are related to certain channels. 

Each set of calibration data comprises the location of the respective 

15 calibration point and at least one measured signal parameter for each of sev- 
eral channels at that calibration point. Calibration points are points whose 
location and signal parameters are known or measured. The calibration meas- 
urements are typically detemnined by means of fixed and/or mobile calibration 
receivers. Fixed calibration receivers can be attached to buildings, traffic signs, 

20 lampposts and the like. IVIobile calibration receivers can be transported with 
persons or in vehicles. The calibration receivers measure the signal parame- 
ters like the actual receivers do. The measured signal parameters can be 
transferred to the statistical model via wired or wireless transmission {=on-line) 
or by moving a detachable media, such as a memory disk, tape or card (=off- 

25 line). 

Location estimation can take place at the receiver site or at the 
network site. If the location is estimated at the receiver site, the receiver (or an 
attached computer) must have access to the statistical model. With cunrent 
technology, a feasible statistical model can be compressed to a size which is 

30 manageable in a laptop or palmtop computer. The model can be updated 
while the computer is connected to the internet, for example. Alternatively, the 
model can be supplied on a detachable memory, such as a CD-ROM or DVD- 
ROM. In the future, even a mobile phone will have sufficient memory for hold- 
ing the statistical model. The model can be updated by means of a data call 

35 via a fast connection, for example. If the receiver site stores a copy of the sta- 
tistical model, It needs no transmission capability, and the actual receiver can 
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be a broadcast reoeiver, a pager or a dedicated add-on card for a laptop com- 
puter, similar in appearance to current GSM attachment cards for laptops. . 

Alternatively, the receiver may be part of a transceiver, such as a 
mobile phone or a \NiAH or Bluetooth interface attached to a portable or 

5 handheld computer. In this case, the transceiver may send the measurement 
results to the network which forwards the results to a location sender. Depend- 
ing on the type of transceiver, the measurements can be sent in a short mes- 
sage, via a data call or a WAP or WLAN connection, for example. The location 
server can send the transceiver Its location estimate over a similar connection. 

10 Accordinig to a prefened emt}odimeht of the invention, the signal 

parameter measurements (the calibration measurements and/or the receiver's 
cun-ent observations) are quantified to a relatively sniall number of classes, 
such as two to five classes. In other words, the granularity of the measure- 
ments is increased. At first sight, such granularity increase seems to lose in- 

15 fpmiation. For instance, assume that the signal strength of a certain channel at 
a certain location is 34 units on a scale of 0 to 100 (the actual unit is inele- 
vantj. Instead of storing the result of 34 units, we only store the fact that the 
measurement was between 25 and 50, ie a value of 1 on a scale of 0 to 3. It 
would seem that a value of 34 on a scale of 0 - 100 can better predict the sig- 

20 nal strength In the neighbourhood of that location than a value of 1 on a scaje 
of 0 to 3 does. However, in many cases increased granularity results in in- 
creased location accuracy. One reason for this is that on a high-resolution 
scale, there are many values that occur relatively seldom, whereas on a low- 
resolution scale, all possible values occur relatively frequently. 

25 An advantage of the invention is that prior infonmation on the net- 

wori< infrastructure is not necessary (although it may be useful). This means 
that a location service according to the invention Is not tied to networic opera- 
tors. Even if the location service according to the invention is maintained by a 
networi< operator, that operator can exploit observations from other operators' 

30 networi^s without prior infonnation on their infrastructure. The invention is ap- 
plicable in a wide variety of networic techniques, such as cellular networi^s, 
broadcast networics or wireless local-area networks. 

Brief description of the drawings 

The invention will be described in more detail by means of prefered 
35 embodiments with reference to the appended drawing wherein: 
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Figure 1 illustrates various graphis of signal parametisr versus re^ 
ceiver location; 

Figure 2 is a block diagram illustrating the general concept of the 
invention; . 

5 Figure 3 is a block diagrarnJIIustFating a typical calibration receiver 

for detennining calibration measurements; 

Figures 4A and 4B arei block diagrams illustrating mobile receivers 
whose location is' to be estimated; and 

Figure 5 illustrates the structure of a statistical model. 

10 Detailed description of the invention 

Figure 1 illustrates various graphs of signal parameter versus re- 
ceiver location. The horizontal axis represents the (one<lirhensionai) location 
of a receiver. The vertical axis represents a signal parameter V (such as signal 
strength, or enror ratio) measured by a receiver. Graphs A and B depict signal 

15 parameters of two channels. In this hypothetical example, we have 10 data 
points Di to Dio measured at location Xi to Xio, respectively. Both graphs A 
and B share the data points Di to Dio having the respective locations Xi to Xio 
and the signal parameter value Vq. Figure 1 gives a feint idea of the difficuities 
in implementing the invention. Not only is the parameter value Vo common to 

20 1 0 different locations (in this example), but the 1 0 locations could be explained 
equally well by both graphs A and B. The well-known Nyquist criterion states 
that a signal can be fully reconstmcted if sampled at more than twice its high- 
est frequency component, if the graphs A and B represent, say. field strength 
in a GSM network having a nominal frequency of 900 MHz, the spatial fre- 

25 quency of the graphs A and B has a wavelength of approximately 30 cm. Ac- 
cordingly, the signal parameters should be sampled at points less than 15 cm 
apart, which is clearly an impossible task. But if the signal parariieters are 
sampled at points more than half a wavelength apart, the graphs A and B can- 
not be reconstmcted. as evidenced by the fact that between points Xe and Xio 

30 the graphs A and B have no similarity whatsoever. 

The reason that the present invention works in practice stems from 
the fact that as the number of channels increases, the number of locations 
where the channels behave as described above decreases rapidly, and so it 
becomes increasingly unlikely that any two points cannot be distinguished 

35 from each other based on the measured parameters. 
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Figure 2 is a blocl^ diagram illustrating the general concept of the 
Invention. In Figure 2, the invention is implemented as a compact location es- 
timation module LEM, although more distributed Implementations are eiqually 
possible. An essential feature of the invention is a statistical modisl SM of the 

5 receiver's wireless environment, tiie model being able to predict tiie receiver's 
location given a plurality of current observations at the receiver site. The statis- 
tical model SM is built and maintained by a model construction module MCM, 
on the basis of calibration data CD anid, optionally, on the basis of prior infor- 
mation Pi of the wireless environment. The optional prior infomnation PI may 

10 comprise Infonnation on networi^ infrastnicture, such as the locations and ra- 
dio parameters of base stations. The locations at which calibration hieasure- 
ments are collected are called calibration points. The calibration data CD 
comprises data records each of which comprises the location X of the calibra- 
tion point in question and the set of signal parameters V measured at that call- - 

15 bration point. Optionally, the calibration data records may also comprise the 
time at which the measurement was made, in case the signal parameters vary 
with time. The location X can be expressed in any absolute or relative coordi- 
nate system: In special cases, such as trains, highways, tunnels, waterways or 
the lii<e, a single coordinate may be sufficient, but nomnally two or three co- 

20 ordlnates will be used. The reference sign X denotes the set of all coordinates 
of the location. 

It should be noted that the tenn 'training data' is often used in the 
context of such statistical models. In tiie context of this invention, the tenn 
'calibration' is prefenred, because 'training' may convey the Idea tiiat the model 

25 is ready after initial training, whereas 'calibration' better conveys the idea that 
the model may have to be updated regulariy as the conditions change. 

There is also a location calculation module LCM for producing a lo- 
cation estimate LE on the basis of the receiver's current observations CO and 
the statistical model SM. Technically, the 'measurements' and 'observations' 

30 can be performed simllariy, but to avoid confusion, the term 'measurement* is 
generally used for the calibration measurements, and the signal parameters 
obtained at the cun-ent location of the receiver are called 'observations'. The 
receiver's most recent set of observations is called current observations. The 
location calculation module LCM or a separate estimate interpretation module 

35 EIM may also use the receiver's observation history OH to interpret the loca- 
tion estimate. In other words, the obsen/ation history OH can be used to re- 
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solve ambiguities in cas^s where a. set of observations can be explained by 
two or (inore locations with substantially equal probability. 

Figure 3 is a block diagram illustrating a typical calibration receiver 
CR for determining the calibration measurements in the calibration data CD 

5 shown in Figure 2. Figure 3 shows a mobile calibration receiver comprising a 
portable computer (or data processor) PC-C, a mobile station MS-C (such as a 
GSI\/I, GPRS or UMTS mobile phone) and a location receiver, such as a GPS 
. (globar positioning system) device. The suffixes -C stand for calibration re- 
ceiver, to distinguish the corresponding parts of the actual receiver R in Figure 

10 4. For clarity, the calibration receiver's main modules PC-C, MS-C and LR are 
shown separately, although the two latter modules are available as PC cards 
which can be inserted into a card socket in a typical laptop computer. The 
calibration receiver CR observes the radio signal parameters of the available 
base stations BS in a cellular radio network RN. The Interface between the ra- 

15 dio network RN and the mobile station MS-C is called a radio interface Rl. If 
the radio interface Rl is bidirectional, the calibration receiver CR may send its 
obsen/ations to the location estimation module LEM via the same radio inter- 
face Rl. Altematively, the calibration receiver's portable computer PC-C may 
store the observations on a detachable memory DM medium, such as a re- 

20 cordable CD-ROM disk, which is later brought off-line to the location estima- 
tion module LEM. 

The location receiver LR of the calibration receiver CR can be en- 
tirely conventional, for example a commercial GPS (global positioning system) 
receiver, as long as it can output the measured coordinates to an attached 

25 computer or other data processor. The portable computer can also be a con- 
ventional, suitably programmed computer. Only the mobile station MS-C may 
heed modifications to its hardware or finnware (its ROM contents). Modifica- 
tions may be needed, depending on how many signal parameters the mobile 
station measures. For example, a conventional GSM phone monitors, in addi- 

30 tion to its currently active cell, some parameters of its neighbouring cells, but 
the neighbouring cells are not measured as extensively as the active cell. Only 
when a GSM phone Is having an active call, does it monitor the neighbouring 
cells as extensively as its active cell. For the purposes of the invention, it 
would be beneficial to modify the mobile station's cell monitoring routines such 

35 that it monitors the available cells as extensively as possible. 
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Naturally, the calibration receiver CR can comprise more than one 
. mobile station for monitoring different types of netwbrlcs or different operator's 
network. For monitoring broadcast networks, the calibration receiver CR 
should aiso comprise a scanning broadcast receiver (not shown separately). 

5 Altematively, the mobile station MS can be a mujti-mode device capable of re- 
ceiving cellular networks and broadcast networks. 

Calibration receivers, like the one shown in Figure 3, can be canied 
along in vehicles or with persons. Fixed calibration receivers, which do not 
need a OP'S receiver, can be attached to buildings, traffic signs, lampposts 

10 and the like. As an altemative to using a separate location receiver, the loca- 
tion of the calibration receiver can be detennined by one or more of the follow- 
ing techniques: showing the receiver's location on a digitized map; entering a 
street (or other) address and converting it to a location by means of a suitable 
database; or using other known locations, such as stops of public vehicles. 

15 Figure 4A is a block diagram illustrating a typical mobile receiver 

whose location is to be estimated. A simple embodiment of a receiver R conv 
prises only a suitably programmed mobile station MS. For some embodiments, 
the receiver R may also comprise a portable computer (or data processor) PC. 
Again, the term 'receiver' implies that the device is receiving when its location 

20 is being estimated although, in practice, most embodiments will also have 
transmitting capability. The embodiment shown in Figure 4A does not contain 
the statistical model SM. Accordingly, the receiver R must send its cunrent ob- 
servation set CO to the location estimation module LEM via the base station 
BS it is connected to. The location estimation module LEM returns the receiver 

25 its location estimate LE via the radio interface Rl. 

Figure 4B shows an altemative embodiment in which the receiver's 
attached computer PC receives a copy of the statistical model SM on a de- 
tachable memory DM, such as a CD-ROM disk, and the receiver is able to de- 
temnine its own location without transmitting anything. As a yet further altema- 

30 tive (not shown separately), the receiver's attached computer PC may receive 
the statistical model via an Internet (or any other data) connection to the loca- 
tion estimation module LEM. Future wideband mobile stations may be able to 
receive the statistical model via the radio interface Rl. A hybrid of the tech- 
nologies may also be used such that the receiver receives an initial statistical 

35 model via a wired connection or on the detachable memory, but later updates 
to the model are sent via the radio interface. 
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Note that in Rgures 3, 4A arid 4B, the radio networic RN Is shown 
as a cellular network and the mobile stations MS resemble Cellular handsets. 
The invention is not limited to cellular networlcs, however, arid can equally well 
be used In a WLAN environment, In which case the mobile statioris ard re- 
5 placed by Wi-AN Interface devices. - 

Statistical modelling 

Possible statistical models will now be studied in more detail. In 
general, a statistical model, as used in this context, can comprise several indi- 
vidual statistical submodels, in which case the actual estimate is obtained by 

10 combining the individual results of the submodels. 

There are many possible statistical modelling approaches that can 
be used for producing the required statistical submodels. In the following we 
will focus on the probabilistic approach. A probabilistic model means that when 
estimating the location of the mobile temninal, the result is represented as a 

15 probability distribution over the possible locations if location X is modelled as a 
discrete variable, whereas, if the location X is modelled as a continuous vari- 
able, the result is represented as a density function, in the following, the focus 
will be on the discrete case. Similariy, the location-dependent measurements 
V can also be modelled either with discrete or continuous observational vari- 

20 ables. The number of dimensions of vector V (the number of measurements 
that can be obtained) varies and depends on the properties of the operating 
wireless network(s). 

Again, there are many probabilistic model classes that can be used. 
In the following prefen^ed embodiment of the invention, the focus will be on pa- 

25 rametric probabilistic models. In this case a single model can be represented 
as a pair (M.O), where M denotes tiie model stmcture, ie the qualitative proper- 
ties of the model that detennine which parameters are required, and 0 denoteis 
the quantitative values of the parameters. 

There are two principal approaches for constructing parametric 

30 probabilistic models (M, 6) in the present context, namely conditional models 
and joint models. Conditional models are models that directly give probability 
distributions of the forrri P(X | V,M, 0), where V denotes the values of the ob- 
servational variables (for example, a vector consisting of signal strength 
measurements), and X denotes the location where obsen/ation V was made. 

35 Joint models define probability distributions P(X,V | iVI, 9) on events (X,V). 
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However, by using the axioms of probability theory we can see that 
P(X I V.IVI, e) = P(X ,V I M, eyP(y | M, e)/where P(V | M, e) does not depend 
on the location X. Thus we can treat the denominator P(V | M, &) as a nonnal- 
izing constant This means that we can always use a joint model for condl- 
5 tional modelling, and in the following we will focus on joint modelling and re- 
gard conditional modelling as a special case. 

There are many ways to use parametric models in location estima- 
tion. Let us first assume that we have decided to use a single model structure 
M, and we wish to detemilne the parameters frbrii the cailbration data CD so 
10 that we get a joint probabilistic model for events (X.V), which also gives, as 
described above, the required conditional distribution for location X given the 
observations V. As described in [Kontkanen et al. 2000], there are several al- 
tematives for producing the joint distribution: 

1. We can use P(X, V | M, e(D)), where e(D) is the maximum likelihood instan- 
15 tiation of the parameters, ie e(D) = arg max P(D 1 1\4, 9). 

2. We can use P(X,V | M, 9(0)), where 9(D) is the Bayesian maximum poste- 
rior instantiation of the parameters, ie 9(D) = arg max P(9 1 M, D) 

3. We can use Pp(,V | M, 9(D)), where 9(D) is the mean of the posterior distri- 
bution P(9 | M. D). 

20 4. We can integrate over the parameters 9: P(X, V | D,M) = fP(X,V | D,M, 9)P(9 
|D,M)d9. 

5. We can use P(X, V | M, 9(D)), where 9(D) is the parameter instantiation op- 
timizing the minimum message length criterion described in [Wallace and 
Dowe 1999] and the references therein. 
25 in some special cases, altematives 3 and 4 are equivalent. 

In general, one may wish to use several model structures M. In the 
following, we will assume that we have fixed the general model family (set) Fi 
the set of all the possible model structures under consideration. For example, 
the set F may conrespond to the set of all possible Bayesian network models 
30 (see [Cowell et al. 1999], [Peari 1988]). In this case we produce the predictive 
distribution P(X |V,F) by computing a weighted sum over all the models in F: 
P(X |V,F) oc S P(X , V I IVI) W(M). Possible weighting functions W include the 
following: 

1. The posterior of the model structure M, given the data: 
35 P(M I D) oc P(D I M)P(M) = P(M) I P(D |9,M)P(9 | M)d9. 
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2, The stochastic complexity of the data, given the mcicle! structure M, and the 
approximations of the stochastic complexity criterion, discussed in 
[Rissanen 1999] and the references therein. 

3. The minimum message length of the data, given the model structure M, and 
5 the approximations of the MML criterion, discussed In [Wallace and Dowe 

1999] and the.references therein. 

It is also possible to use conditional (supervised) versions of the 
weighting functions, in which case the weights are computed with respect to 
conditionar modelling, and the actual data is taken to consist of only the values 

10 of the location variable X, and the measurement data V is treated as "back- 
ground data", these alternatives are discussed in [Kontkanen et al 1999J. 

If the number of model structures M in F is too high for computing 
the weighted sum In a feasible time, we have to restrict the model family F by 
perfonnlng a search in F, and pruning F to consist of only those model stmo- 

15 tures that are the best with respect to some cost function. The possible cost 
functions for perfomiing the search include the weight functions listed above. 
Any search algorithm can be exploited In this task. An extreme case of this 
type of restricting search is a case Where only one single model stmcture M in 
F is chosen, in other words, the sum over model structures reduces to a single 

20 term conresponding to the use of a single model with the largest weight. 

If the observations V are modelled as discrete variables, the granu- 
larity of the discrete variables can be viewed as part of the model structure M. 
The granularity can either be fixed by the user (representing prior Infonmation), 
or as pari of the model structure M, It can be learned from the calibration data. 

25 The optional prior infomiation, such as information on the locations 

and radio parameters of the base stations, represents knowledge other than 
that extracted from the calibration measurements. In the probabilistic setting, 
we can Identify the following ways for coding the prior Information: 

1. By choosing the initial model family F of probability models (detennlning the 
30 model structures considered, and with each model structure, the fomns of 

the distributions used and the assumptions made). 

2. If the observational variables V are taken to be discrete, by choosing the 
granularity of the discretization. 

3. If the location variable X is taken to be discrete, by choosing the granularity 
35 of the discretization. 
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4. By determining the prior distribution P(e | M) for the parameters of the 
model M. 

is. By determining the prior distribution P(M) for the model structurfes M In the 
family F. 

5 Missing data 

There are several altematlve procedures for handling missing data: 

1 . Treat 'missing' as an extra value for the variable In question. 

2. Ignore the missing entries (the sufficient statistics are computed from the 
existing data only) 

10 3. Estimate the missing values from the existing data and/or prior infonnation. 
The estimates can either be used for filling in educated guesses of the 
missing values, or they can be treated as partial observations (sufficient sta- 
tistics of several possible values can be simultaneously partially updated, 
according to, e.g., their estimated probabilities). 

15 4. Fill In the missing values by using random guesses. 

Location interpretation and reporting 

The result of the probabilistic location estimation can be reported in 
several different ways. First, we can divide the working area Into several sub- 
areas in different ways: the subareas can either form a full partitioning of the 
20 wori< area, or they can cover only a portion of the whole work area. An exam- 
ple of the latter case is that only the locations listed in the calibration data D 
(with a desired accuracy) are considered. The result of the probabilistic loca- 
tion estimation can now be reported either 

1 . By giving the full probability distribution over the areas, ie, for each area X, 
25 give the conresponding probability P(X I V,F). 

2. By giving the most probable subarea X with respect to the distribution P(X | 
V.F). 

3. By giving a point estimate minimizing the expected value of some error 
function with respect to the distribution P(X | V,F). 

30 An example of altemative 3 Is the mean squared enor, In which 

case the point estimate Is the vveighted average of the centre points of the su- 
bareas (assuming that the subareas are of equal size), the weights being the 
probabilities P(X | V,F). If the subareas X are not of equal size, the weights 
can be rescaled with respect to the relative size of the conresponding subarea, 

35 for example by multiplication. 
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Uncertainty about the receiver's location can be reduced by prior in- 
fomiation PI, if available, and/or the observation history OH. Let us assume 
that the above alternative 1 was chosen Initially. In other words, the user or 
application requesting the location of the receiver is reported the full probability 

5 distribution. The probability distribution may indicate a number of feasible loca- 
tions. The prior Infonnation PI, if available, may indicate that only one of the 
locations is possible, given the received cell identifiers or the like. Alternatively, 
the observation history OH can be used to exclude some locations. For exam- 
ple, although a number of locations could explain the receiver's current loca- 

10 tion, only a subset of the locations can explain the entire observation history 
OH, given the receiver's finite speed. 

Performance examples 

Example 1 : Location estimation with the NalVe Bayes model. 

The subareas X under consideration are the locations where the 

15 calibration data was collected. The radius of the locations is assumed to be 
one meter, although any unit can be used. The observational variables V are 
taken to be discrete with m values. The value of m can be a constant (eg 3), or 
it can be optimized by using one of the weighting functions described above. 
The boundary points between the intervals can be determined so that the 

20 number of training samples within each interval is the same (equal-frequency 
discretization), or alternatively, the Intervals can be made of equal width 
(equal-width discretization). The intervals can also be detemiined by using a 
clustering algorithm, such as the K-means algorithm. 

One model stmcture M is used: the observational variables Vi,...,Vn 

25 are . assumed to be independent, given the value of the location variable X. 
The data is assumed to be independently and identically distributed (= "i.i.d."), 
and follow the Multinomial distribution with the assumptions described in [Gei- 
ger and Heckemnan, 1998], Prior information is non-existent A non- 
informative uniform prior distribution for the model parameters is used. 

30 Alternatives are discussed in [Kontkanen et al, 2000]. The distribution P(X,V | 
D,M) is computed by integrating over the parameters. With the above 
assumptions, this can be done as described in [Kontkanen et al., 2000]. 

In this experiment, the obsen^ation history OH is taken into account 
by treating the eight (other numbers are equally possible) last signal meas- 

35 urements as a isingle measurement vector V so that the eight individual meas- 
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uremehts are assumed to be independent of each other. The result is given as 
a point computed as a weighted average of the centre points of the subareas, 

where the weight for area X is P(X I V.b.M): 

This method was implemented and tested empiricaiiy in downtown 

5 IHelsinl(i by exploiting the signal strengths of Sonera GSM network. The work 
area was approximately 400 x 500 meters in size, and the calibration data was 
collected outside |n the streets in approximately 50 evenly distributed points. 
The average distance, between two measurement locations was approximately 
. 50 nieiers. The system was tested by using the location estimator in 300 ran- 

10 domly situated locations within the work area. The average location enor in 
this test was 42 meters. 

Example 2: Location estimation with a Mixture of Histograms model. 

The location variable X is taken to consist of two coordinates (one 
or three coordinates are also possible) on a fine-grained, discrete scale. The 

15 resolution of the scale can be, for instance, one meter. The observational vari- 
ables Vi, Vn are taken to be discrete with the maximum resolution deter- 
mined by the measuring device, e.g. 1 dBm. We denote the combination of Vi, 
Vn by V. Missing values are replaced by a value which is smaller than any 
possible observable value. Several models are considered. Each model Mm is 

20 associated with parameters k, I, and Gid whose semantics is described below. 

Figure 5 illustrates the stmcture of model Mm. The value of variable 
Xk is obtained from the value of variable X via discretization into k values. The 
conditional distribution of variables Vi(l), Vn(l) given the value of Xk is de- 
scribed by model parameters Gw . Each Vi, where I belongs to the set {1 , .... n}, 

25 is unifomily distributed within the interval defined by the value of variable Vi(l). 
A low-resolution location variable Xk Is derived from the fine-grained location 
variable X by discretizing the latter with equal width discretization (other discre- 
ti^tion methods also possible) using k bins, le k possible values. Whenever a 
boundary point of the low-resolution discretization appears between two 

30 boundary points of the fine-grained discretization, the mass (i.e. the number of 
observations within the sub-interval) is divided according to the relative size of 
the overiapping intervals. For instance, let the fine-grained discretization have 
5 bins (4 boundary points) within the range [0, 10]. Let the low-resolution dis- 
cretization have 2 bins, and therefore one boundary point at the value 5. If 

35 there are n obsen/ations within the range [4, 6], i.e. within the third bin of the 
fine-grained discretization, then both low-resolution bins get n/2 observations, 
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because the boundary point 5 split^ the range [4, 6] into two parts of equal 
size. Lilcewise, each observational variable. V| is discretized using I possible 
values, thus obtaining a low-resolution variable V|(I). 

Model Mki describes the conditional probability functions P(V(I) | Xk. 

5 Mm, Old), where Oid denotes the model parameters of model Mm. The low- 
resolution observational variables Vi(l), Vn(l) are taken to be independent, 
given the value of the location variable Xk. For each i belonging to the set {1 , 
n}, the distribution P(Vi{l) | Xk, Mm, Om) Is talcen to be i.i.d., and follow ttie 
Muitinomiai-Dirichjet distributibri witti the assumptions described in [Geiger 

10 and Heckennan, 1998]. Prior infonmation Is nonexistent. A unifonnn prior distri- 
bution over models Mid is used. A non-infonmative equivalent sample sjze 
(ESS) prior distribution for the model parameters is used as described in 
[Heckennan, 1995]. A second-order prior for the ESS parameter Is used, e.g., 
a uniform distribution over tiie set {1 , 1 0}. For isach 1 belonging to tfie set {1 , . . . , 

15 h}, the distribution P(V|(I) | Xk, Mm) is computed by integrating oyer the model 
parameters. With the above assumptions, this can be done as described in 
[Konti^anen et al., 2000]. 

The distiibution P(V| { Vi(l)) is taken to be uniform over the interval 
defined by the value of Vi(l) and the discretization of V| defined by parameter I. 

20 For instance, let the range of V| be [0, 10], let ttie value of I be 5, and let the 
value of Vi(l) be 2. Assuming that equal widtti discretization is used, the values 
of Vi are discretized into five intervals, [0, 2], [2, 4], [4, 6], [6, 8], and [8, 10]. 
Now, given that the value of Vi(l) is 2, the distribution P(V| | Vi(l)) is unifomrt 
over the interval [2, 4]. The variables Vi, Vn are taken to be independent of 

25 each other given the Values of variables Vi(l), Vn(l). 

Combining the two distributions P(V(I) | Xk, Mm) and P(V | V(l)), we 
obtain a conditional distribution P(V | Xk. Mm). The distribution P(V | X, D) is 
computed as a weighted average over the models Mm, where k and I vary over 
ttie set {1, 20} (other choices equally possible). The models are weighted 

30 by the marginal likelihood P(V(D) | Xk(D), Mm), where the calibration data (with 
n observations) consists of the vectors V(D) = (V^D), .... V"(D)) and Xk(D) = 

(Xk'(D) Xk''(D)). 

With these assumptions, the marginal likelihood can be computed 
efficiently in two parts: First, the product of tiie tenns with the fomi P(V(I) | Xk. 

35 Mm) can be computed as described in [Heckennan, 1995] and [Geiger and 
Heckerman, 1998]; Second, the tenns with ttie form P(V | V(l)) have the same 
value, which is a constant depending on I, because the distribution P(V | V(l)) 
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is uniform. The result is givien as a posterior probability distribution P(X | V, D) 
= P(V I X, D). P(X I pj / P{V I i?) over tlie location variable X. The distribution 
P(X I D) is tal(en to be unifonri; Tiie temi P(V | D) is a nomnalizlng factor 
whose value is ignored. Instead the resulting distribution P(X | V, D) is nofmal- 

5 ized so that it sums up to one. 

The method described above was implemented and tested empiri- 
cally in Helsinki on the seicond floor of the building at address Teoliisuuskatu 
23 by using a laptop computer measuring the WI.AN signal strengths through 
a WLAN PC cani. The woric area was approximately 20 x,45 meters (900 

10 square; meters) in size. Calibration data was collected in 12 arbitrary places, 
and the total number of data vectors in it was 204. The system was tested by 
using the location estimator in 25 randomly chosen locations within the woric 
area. Location estimation was repeated five times in each location. When tiie 
above system was used for determining the location area with 95% of the 

15 probability mass, the cpnrect place was in this area 77 % of the time. The av- 
erage size of the 95 % probability mass area was approximately 151 square 
meters, i.e. about 1 7 % of the total area. 
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Claims 

1. A method ifbr estimating a location (X) of a receiver (R, R') in a 
wireless telecommunication environment (RN), the telecommunication envi- 
ronment comprising a plurality of channels for simultaneous communication, 

5 each channel having at least one signal paranieter (V) that varies with location 
(X) differently from the other channels; 

characterized in that the method comprises the steps of: 
for each of a plurality of Calibration points In the wjreless telecom- 
munication environment, detemnining a set of calibration data (CD), each set 
10 of calibration data comprising the location (X) of the respective calibration 
point and at least one measured signal parameter (V) for each of several 
channels at that calibration point; 

on the basis of the sets of calibration data (CD), maintaining a sta- 
tistical model (SM) of the signal parameters (V) of the several channels versus 
15 a receiver's location in the wireless telecommunication environment (RN); 

determining a set of observed signal parameters (CO), the set com- 
prising at least one observed signal parameter (V) for each of several chari- 
nels at the location (X) of the receiver (R, R'); and 

determining a location estimate (LE) approximating the location (X) 
20 of the receiver (R, R') on the basis of the statistical model (SIVI) and the set of 
observed signal parameters (CO). 

2. A method according to claim 1, characterized by the re- 
ceiver (R) sending the set of observed signal parametere (CO) to an external 
location estimation module (LEM) which sends the location estimate (LE) to 

25 the receiver. 

3. A method according to claim 1, characterized by the re- 
ceiver (R') storing a copy of the statistical model (SM) and determining the lo- 
cation estimate (LE) on the basis of the copy of the statistical model (SM). 

4. A method according to any one of the preceding claims, cha- 
30 racterized by maintaining the statistical model (SM) also on the basis of 

prior information (Pi) about the wireless environment's (RN) infrastructure. 

5. A method according to any one of the claims 1 to 4, characte- 
rized in that the statistical model (SM) Is or comprisies a probabilistic model, 
preferably a Bayesian model. 
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6. A method according to clairn 5, c h a ra ct e r ize d in that the sta- 
tistical model (SM) is or comprises a Bayesiati network model. 

7. A method according to any one of the preceding daims, cha- 
racterized in that the signal parameters (V) in the statistical model (SM) 

5 are independent of each other, given the location (X). 

8. A method according to any one of the preceding claims, cha- 
re cteri zed by reducing uncertainty concerning the receiver's location on the 
basis of a history (OH) of the observed signal parameters. 

9. A method according to any of the preceding claims, c h a r a c t e - 
10 rized by rtiodelling at least some of the signal parameters (V) by discrete 

variables whose values conespiond to intervals or unions of intervals on the 
range of possible signal parameter values. 

10. A method according to any of the preceding claims, charac- 
terized by modelling the location (X) as a discrete variable. 

15 11. A location estimating apparatus (IBM) for estimating a location 

(X) of a receiver (R, R') in a wireless telecommunication environment (RN), the 
telecommunication environment comprising a plurality of channels for simulta- 
neous communication, each channel having at least one signal parameter (V) 
that varies with location (X) differently from the other channels; 

20 characterized by: 

a model construction module (MCM) for 

- receiving a set of calibration data (CD) for each of a plurality of 
calibration points in the wireless telecommunication environment, each set of 
calibration data comprising the location (X) of the respective calibration point 

25 and at least one measured signal parameter (V) for each of several channels 
at that calibration point; and 

- maintaining, on the basis of the sets of calibration data (CD), a 
statistical model (SM) of the signal parameters (V) of the several channels 
versus a receiver's location in the wireless telecommunication environment 

30 (RN); 

and a location calculation module (LCM) for: 

- receiving a set of observed signal parameters (CO), the set com- 
prising at least one observed signal parameter (V) for each of several chan- 
nels at the location (X) of the receiver (R, R'); and 
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- dbtermihing a location estimate (LE). approximating ttiis location 
(X) of tlie receiver (R, R') on the basis of the statistical model (SM) and the set 
of observed signal parameters (CO). 

12. A receiver (R, R') comprising means for detemnining sets of ob- 
5 served signal parameters (CO), each set comprising at least one obsen/ed 

signal parameter (V) for each of several channels at the location (X) of the re- 
ceiver (R), characterized by means for conveying the sets of observed 
signal parameters (CO) to a location calculation module (LCM) for detennining 
a location estimate (LE) approximating the location (X) of the receiver (R) on 
10 the basis of said sets and a statistical model (SiVi) of the signal parameters (V) 
of the several channels versus a receiver's location in a wireless telecommuni- 
cation environment (RN). 

13. A receiver (R') according to claim 12, characterized by 
comprising the location calculation module (LCIVI). 

15 14. A receiver (R) according to claim 12, characterized in that 

the means for conveying the sets of observed signal parameters conriprises 
means (Rl) for conveying the sets to an extemal location calculation module 
{LCM). 

15. A receiver according to any one of claim 12 to 14, characte- 
20 rized in that at least some of the sets of observed signal parameters (CO) 
relate to networics the receiver is not attached to. 



wo 02/054813 



I Signal parameter (V) 




xe X7 X8 



Calibration 
measurements X, V, (T) 






CD CO 




Statistical model 



Calculation rules 
and parameters 



CR 







> 


r 


Location 
calcuation (X) 



Location 



LCM 



-1 1 



estimate ! Estimate n . 
r~^j interpretation rT** 
) — y— 'I 



LE 



LEM 



EIM 




wo 02/054813 



PCT/nOl/OlliSl 




wo 02/054813 



INTERNATIONAL SEARCH REPORT 



PCT/FI 01/01151 



A. CIASSIFICATIGN OF SUBJECT MATTER 

IPC7: H04Q 7/38, 601S 5/02 . ' 

According to International Patent Classification (IPC) or to both nflUonal dassification and IPC 



B. FIELDS SEARCHED 



Mit±num documentation searched (classincation system followed by classincation symbols) 
IPC7: H04Q 



Documentation searched other than minimum documentation to the extent that such documents are included in the fields searched 

SE,DK,FI,NO classes as above 



Electronic data base consulted dudng the international search (name of data base andi where.practicable, seardi terms used) 



C. DOCUMENTS CONSIDERED TO BE RELEVANT 



Category* 



Citation of document, with indication, where appropriate, of the relevant passages 



Relevant to claim No. 



US 6112095 A (WAX ET AL), 29 August 2000 

(29. 08 • 00), column 7, line 24 - column 8/ line 42, 
abstract 



EP 1045602 A2 (LUCENT TECHNOLOGIES INC.), 
18 October 2000 (18.10.00), column 3, 
line 45 - column 4, line 14; column 4, 
line 34 - line 49, summary 



US 5481588 A (RICKLI ET AL), 2 January 1996 
(02.01.96), column 2, line 45 - line 59 



1-15 



1-15 



1-15 



~)(| Further documents are listed m the continuation of Box C. |- xl See patent family annex. 



* Special categories of cited documents: 

"A" document deGning the genera] state of the art which is not considered 

to be of particular relevance 
"E" earlier application or patent but published on or after the interoatianal 

filing date 

"If document whi<fh may doubts on priority claim(s) or «4uGh is 
dted to establish the publication date of anottier citation or other 
special reason (as s]}edfied) 

'O' document referring to an oral disclosure, use, exhibition or other 
means 

'P'' document published prior to the international filing date butlater Uian 
the priority date claimed 



*T' later document publidied after the intemational filing date or priority 
date and not in conflict with the application but dted to understand 
the principle or theory underlying me invention 

'X' document of particular relevance: the claimed invention cannot be 
considered novel or cannot be considered to invdve an inventive 
step when the document is taiken alone 

'Y* document of particular relevance: the claimed invention cannot be 
considered to involve an inventive step when the document is 
combined with one or more other su^ documents, such combination 
being obvious to a person skilled in the art 

document member of the same patent family 



Date of the actual completion of the intemational search 



12 April 2002 



Date of mailing of the intemational search report 

A 1 -04- 2002 



Name and mailing address of the ISA/ 
Swedish Patent Office 
Box 5055, S-102 42 STOCKHOLM 



Authorized officer 



Mai In Soderman/js 

4 Xt. .1 MK n p^m 



wo 02/054813 

INTERNATIONAL SEARCH REPORT 


International .a{^ucauon nu. 
PCT/FI 01/01151 


C (ConUnuaeon); DOCUMENTS CONSIDIREP TO BE RBE^E^ 


Category* 


GtaUon of* docuineat,.«ijlh indication,. where appropriate, of tlie releraat passages 


Retevant to daiin No. 


P,A 


US 6263208 Bl (CHANG ET AL), 17 July 2001 
(17.07.01), claim 1, abstract 




1-15 


P.A 


US 6266014 Bl (FAHOUCHE ET AL), 24 July 2001 
(24.07.01), claims 1,2, abstract 


1-15 



WO02/y54813. 

information on patent family members 28/01/02 



International ap^,^^^,j,pjj5j 

PCT/FI 01/01151 



Patent docniment 

. dted in search report 


Pub]icaUon 
date - 


pittent family 
]^ember(s) 


Publication 
date 


us 6lir2095 


A 




29/08/00 


US 


6026304 A ; 








us 


6(i64t39 A 










us 


6084546 A 












US 


6101390 A 


08/08/00 










us 


6108557 A 


22/08/00 










us 


6232918 B 


ic/nc/ni 










us 


6249680 B 


i!7/Up/UX 


EP 1045602 


A2 


18/10/00 


AU 


2641000 A 


iQ/in/nn 






BR 


0004836 A 


^U/Uc/UJL 










CA 


2303633 A 


14/10/00 










CN 


1270462 A 


18/10/00 










JP 


2000324537 A 




• ■ ■ • • 








US 


633042d B 


11/10 /ftl 


US 5481588 


A 




02/01/96 


AT 


157835 T 


1 C /AO /07 










AU 


674384 B 


1 Q /I 0 /QC 










AU 


5920594 A 


1 o /I n /Q^ 










CA 


2120580 A 


rt7 /I ri /dii 

u//iu/y4 










DE 


59403915 D 


AA /AA /AA 

yu/uu/uu 












619687 T 


1 O /I A /07 

lo/iu/y/ 










EP 


0619687 A,B 


12/10/94 










SE 


0619687 T3 












ES 


2108931 T 


01/01/98 










NO 


940977 A 


05/10/94 










NZ 


260220 A 


24/04/^7 


US 6263208 


Bl 


17/07/01 


NONE 




US 6266014 


Bl 


24/07/01 


WO 


0165271 A 


07/09/01 



Pnrm PfTr/iRA/7in foatent famiW annex> fJiilv 199X^ 



This Page is Inserted by IFW Indexing and Scanning 
Operations and is not part of the Official Record 



Defective images within this document are accurate representations of the original 
documents submitted by the applicant. 

Defects in the images include but are not limited to the items checked: 

□ BLACK BORDERS 

□ IMAGE CUT OFF AT TOP, BOTTOM OR SIDES 

□ FADED TEXT OR DRAWING 

□ BLURRED OR ILLEGIBLE TEXT OR DRAWING 

□ SKEWED/SLANTED IMAGES 

□ COLOR OR BLACK AND WHITE PHOTOGRAPHS 

□ GRAY SCALE DOCUMENTS 



□ REFERENCE(S) OR EXHIBIT(S) SUBMITTED ARE POOR QUALITY 

□ OTHER: 



IMAGES ARE BEST AVAILABLE COPY. 
As rescanning these documents will not correct the image 
problems checked, please do not report these problems to 
the IFW Image Problem Mailbox. 



BEST AVAILABLE IMAGES 




LINES OR MARKS ON ORIGINAL DOCUMENT 



